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Abstract 

We report a study of semileptonic B decays to P-wave D** mesons. Semileptonic decay to 
a D\ meson is observed for the first time and its product branching ratio is measured to be 
B(B+ -► Dfe+v) x B(Df -»• D-TT+) = 0.22 ± 0.03(stat.) ± 0.04(syst.)%. The result is obtained 
using fully reconstructed B tags from a data sample that contains 657 x W 6 BB pairs collected at 
the T(4S) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. 

PACS numbers: 13.20.-v, 13.20.He, 14.40.Lb 
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Heavy Quark Effective Theory (HQET) has proven to be very successful at describing 
semileptonic decays of P-mesons, especially inclusive transitions; it allows one to extract 
\Vcb\ to better than 2% accuracy However, some difficulties arise when it is applied to 
exclusive decays. For example, certain sum rules (in particular, the Uraltsev sum rule 0) 
imply the strong dominance of decays to the narrow excited D mesons over those to the wide 
excited D mesons, while some experimental data show the opposite trend 0, H| . However, 
no complete experimental study of such semileptonic decays to excited P-mesons exists, and 
thus no direct comparison with theoretical predictions can be performed. In this paper we 
report on a study of B — > D^nlv decays and measure the excited D contributions to the 
D^tt final state. 

According to HQET there are two doublets of orbitally excited (P-wave) charmed mesons 
(D**), differentiated by their light quark angular momentum j q = 1/2 or j q = 3/2. Members 
of the j q = 3/2 doublet are predicted to decay only via a P-wave and be relatively narrow, 
while members of the j q = 1/2 doublet are predicted to decay only via an S*-wave and be 
relatively broad [5]. The P** states with spin-parity and light quark angular momentum 
combinations + (j q = 1/2), l+(j q = 1/2), l+(j q = 3/2) and 2+(j q = 3/2) are usually 
labelled Dq, D[, D\ and D\, respectively. The D** states have previously been observed 
and studied in hadronic P-decays Semileptonic P decays to narrow D\ and D\ mesons 
have been studied by a number of experiments [3]. The semileptonic branching fractions of 
B — > D^islv decays were recently measured by Belle [s[ and BaBar 

This measurement is based on a data sample that contains 657 million BB pairs, which 
corresponds to 605 fb -1 , collected at the T(45') resonance with the Belle detector operating 
at the KEKB asymmetric-energy e + e~ collider [Ic| . An additional 68 fb -1 data sample taken 
at a center-of-mass energy 60MeV below the T(45') resonance is used to study continuum 
e + e~ — > qq (q = u, d, s, c) background. The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement 
of time-of-ffight scintillation counters (TOF), and an electromagnetic calorimeter comprised 
of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that provides a 
1.5 T magnetic field. An iron flux- return located outside the coil is instrumented to detect 
K\ mesons and to identify muons (KLM). The detector is described in detail elsewhere (llj . 
Two inner detector configurations were used. A 2.0 cm beam pipe and a 3-layer silicon 
vertex detector was used for the first sample of 152 million BB pairs, while a 1.5 cm beam 
pipe, a 4-layer silicon detector and a small-cell inner drift chamber were used to record the 



remaining 504 million BB pairs [12 



To suppress the large combinatorial background expected in the reconstruction of final 
states including a neutrino, we use a full reconstruction tagging method, which has been 
improved in comparison with our previous paper 0. The first B meson (denoted as B s {) 
is reconstructed in the semileptonic mode of interest, i.e. as a combination of all final 
particles D^ni except for the neutrino. The remainder of the event is combined into either 
a 7r =l= (n < 6) or D^*' p~ combination to form the tagging B meson (referred to below as 
-Stag)- Semileptonic decays are identified by a peak around zero in the missing mass squared 
spectrum, = (Pbeams — -Pag — -Pi) 2 , where Pbeams is the total four-momentum of the beams 
and P tag and P s \ are the reconstructed four-momenta of the P s i and B tag , respectively. This 
method provides significantly improved resolution in the missing momentum in comparison 
with non-tagging methods, thus allowing background suppression, separation of different 
decay modes and precise calculation of the decay kinematics. 
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Charged tracks are required to originate from the interaction point (IP). Charged tracks 
positively identified as electrons or muons with \P\ > 1.0 GeV/c are used as leptons. Charged 
kaons are identified by combining information on track ionization loss (dE/dx), Cherenkov 
light yields and time-of-flight information. No pion identification is required. Photons are 
identified as isolated electromagnetic showers with energies greater than 50 MeV that are 
not matched to any charged track. 

Kg mesons are reconstructed from 7r + 7r~ pairs having an invariant mass within 
±30MeV/c 2 (±5cr) of the nominal K s mass and a well reconstructed vertex displaced 
from the IP. tt mesons are reconstructed from 7 pairs having an invariant mass within 
±15MeV/c 2 (±3cr) of the tt° nominal mass. Such combinations are then fitted with a tc 
mass constraint to improve the momentum resolution. p + mesons are reconstructed from 
7r + 7T° pairs having an invariant mass within ±0.3GeV/c 2 of the nominal p + mass. 

We reconstruct D° candidates using six decay modes (l3| : K~tt + , K^tt + tt , K~tt + 7t + tt~ , 
K s 7t + tt~, K~K + and K s 7i°, and D + candidates using decays to K^it + , K s tt + 7t + tt~ , 
K~tt + 'k + and K + K~tt + . A ±15MeV/c 2 interval around the appropriate nominal D mass 
is used for all modes except for D° — > K~7r + 7r°, where we use a ±25MeV/c 2 window, 
corresponding to about ±3cr in each case. Selected D candidates are then subjected to a 
mass-vertex constrained fit to improve their momentum resolution. D*° candidates are re- 
constructed via decays to D°7r° and -D°7; the latter mode is not used for B s \ reconstruction 
because in this case the D*tc invariant mass spectrum has a large background. D* + candi- 
dates are reconstructed in two modes: D°ir + and D + tt°. The mass difference M(D*) — M(D) 
for D* candidates is required to be within ±3MeV/c 2 (for Dn) and ±10MeV/c 2 (for D 0r y) 
intervals around the mo* — itid nominal value (about 3cr and 2a, respectively). Recon- 
structed D* candidates are subjected to a mass-vertex constrained fit. 

For the signal B s \ meson side we form Di (normalization mode) and D^iri (signal 
modes) combinations. The energy difference AE = E t3ig — Eqm and beam-constrained 

mass M bc = \J Eq M — P t 2 g , (where E ta _ g and P tag are the tag B candidate center-of-mass 
(CM) energy and momentum and E CM = y/s/2 ~ 5.29 GeV) variables are used for B tag 
selection. The B tllg signal region is defined as M bc > 5.27GeV/c 2 , \AE\ < 40 MeV, which is 
about 3d in both cases. The B tlLg candidates are subjected to an energy constrained fit to 
improve P tag resolution. We also use events from the 50 MeV < \AE\ < 130 MeV sidebands 
for background subtraction. In these events the B tag candidate is fitted with its energy con- 
strained to the center of the sideband. In the case of multiple entries in the signal region, 
the -Btag B a \ candidate with the minimum xiot * s chosen, where xLt ^ s calculated as a sum 
over the \ 2 °f intermediate D and D* mesons mass (mass-vertex) fits and x 2 °f the B tag 
energy fit. The same single candidate selection is applied to sidebands. The average number 
of candidates per event is 1.3 in each case. 

The M 2 spectra for the four semileptonic decays B — ► D^niu are shown in Figs. [U, 
la)-ld) as points with error bars. Clear peaks are evident in all distributions. 

We divide the backgrounds into the following categories: 

(1) Continuum e + e~ — > qq events. 

(2) Backgrounds with the B tag misreconstructed from particles belonging to the other B 
meson or fake tracks. 

(3) B s \ backgrounds with the _B tag reconstructed correctly, which can be further separated 
by their source: 
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(3a) Combinatorial background under the signal from B s i. 
(3b) Hadrons misidentified as leptons. 

(3c) Feed-down from B — > D*-nlv reconstructed as B — > Di^iv with lost neutral(s). 

All backgrounds except for (3c) are reliably determined and finally subtracted directly 
from the data. Backgrounds (1) and (2) are estimated using AE sidebands. Continuum data 
and generic BB MC simulation show that these backgrounds have flat AE distributions, thus 
justifying this procedure. Background (3a) is subtracted using (M(Z)W)) sidebands. 
However the region where backgrounds (2) and (3a) overlap is subtracted twice by using AE 
and M(D^) sidebands. To account for this over-subtraction we use two-dimensional AE 
and M(D&) sidebands. The sum of Ml distributions from one-dimensional M(D^*>) and 
AE sidebands after subtraction of the two-dimensional sidebands are illustrated in Figs. HJ 
la)-ld) by the hatched histograms. These histograms represent the sum of backgrounds 
(1), (2) and (3a). 

Background (3b) is studied with data by using combinations of D^n with high momenta 
hadrons {h + ), where the h + candidate is selected with a lepton veto requirement. The com- 
binatorial backgrounds are subtracted from the observed distributions using AE and 
M(D^*>) sidebands. The obtained M 2 spectra are then multiplied by the known misiden- 



tification rate, which depends on the hadron laboratory momentum (see |14| for details). 
It is found that the remaining small peak around zero, which is due to the contribution of 
B — > D ( -*^nn + (n°) decays, is quite small (~ 0.5 — 1.0%) and we ignore it below, including it 
CIS db systematic error. 

Background (3c) is observed only in the B — > Di^lv channels and is estimated from a 
MC simulation with normalization fixed to the data using B — ► D*(jr)£v signal yields. This 
contribution is plotted in Figs, dj la), 1c) as open histograms. 

The background-subtracted distributions are shown in Figs, [fl 2a) -2d). These dis- 
tributions are fitted with signal functions, the shapes of which are fixed from MC studies. 
Fitted signal yields, reconstruction efficiencies and branching ratios are summarized in Ta- 
ble HI The branching ratios are calculated relative to the normalization modes B — > Dlv to 
cancel out the B taig reconstruction efficiency according to the formula: 

B(mode) = B(norm) x — x 



^norm ^mode 

where iV norm ( mo de) and e nor m(modc) are the signal yield and reconstruction efficiency of the 
normalization mode (mode of interest) and the normalization mode E is taken from the 



PDG [15j. Relative efficiencies are obtained from MC simulation. Intermediate branching 
fractions are included, while the tagging efficiency is not. The reconstruction and back- 
ground subtraction procedures for the B — > Dtv mode are identical to those applied for the 
studied channels. The obtained branchingfractions are in good agreement with our previous 
measurement and with BaBar results |0] • The low efficiency in the last mode is the result 
of not using the D*° — > D ^ decay channel. 

Signals for semileptonic B decays to orbitally excited D** are extracted from the D^n 
invariant mass distributions. We define a signal window for B — > D^niv decays by the 
requirement \M%\ < 0.1 GeV 2 /c 4 . The backgrounds are estimated in the same way as in the 
Ml distribution study. The D^ic invariant mass spectra from the signal window after sub- 
traction of backgrounds (1-3) are shown in Fig. [21 The mass distributions before background 
subtraction, restricted to the region near the j q = 3/2 states, are shown in the insets. 
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TABLE I: Results for B -> D^wiv where the first error is statistical and the second is systematic. 



Mode 


Yield 


Eff.,% 


fi(mode),% 




2320 ± 60 


6.4 


2.15 ± 0.22 a 


B + -» D--K+1+V 


192 ± 19 


2.8 


0.40 ± 0.04 ± 0.06 


B+ D*--k + 1 + v 


123 ± 14 


1.14 


0.64 ±0.08 ±0.09 


B° -► D~i + v 


760 ± 30 


3.7 


2.12 ±0.20 a 


B° -► D°tt-£ + u 


150 ± 20 


3.7 


0.42 ±0.07 ±0.06 


B° -» D*°ir~e+is 


22 ±8 


0.40 


0.56 ± 0.21 ± 0.08 



°Used as a reference. 



To extract the D** signals we perform simultaneous unbinned likelihood fits to the signal 
and background D^n mass spectra. The signal function includes all orbitally excited D** 
contributing to the given final state (D and D\ to Dn and Di, D[, D\ to D*tt), each of 
which is described by a relativistic Breit-Wigner function for a known orbital momenta, and 



a non-resonant part described by the Goity- Roberts model [16| • D** masses and widths are 
fixed to measured values 0. To further investigate the Dir mass spectrum we also test a 
D* + D\ hypothesis. Despite the D°tt + mass region corresponding to D* + being excluded 
from the study, and while D*° is below the D~ir + threshold, a virtual D* can be produced 
off-shell. We describe the D* contribution by a tail of the Breit-Wigner function with floating 
normalization. Fit results are shown as a dashed line for this combination. 

A study of the sidebands shows that the background is described by the sum of a signal 
function and an exponential. The resulting signal function and contributions from the 
resonances are shown in Fig. [2] as solid and dashed curves, respectively, superimposed on 
the background-subtracted mass spectra. In the insets the solid and dashed curves represent 
the fitted signal and background, respectively. In B° — ► D*irlv decays a small feature may 
be observed around 2.6GeV/c 2 , which is absent in B + — > D*n£v. However, the significance 
of this feature is small (2.5a) and there is no known state there, so we do not include a term 
for it in the fit. Fitted resonance yields and corresponding product branching ratios are 
listed in Table [TTJ The contribution of the non-resonant component in all cases is consistent 



with zero. The B — > D**£u decay significance is defined as y— 2 lnL max /L , where L is the 

likelihood value returned by the fit to the D^n distribution with the D** contribution fixed 
to zero. Our result for B — > Di£ + u is in good agreement with previous measurements Q. 
For a + hypothesis the branching ratio of the decay to the wide Dq is large, in contrast 
to theoretical predictions However, the present statistics do not definitely exclude an 
interpretation of broadly distributed Dtt + events as the D* tail. 

For .D*'** 's decaying into Dtt we perform a study of the helicity angle distributions, which 
is the angle between 7r momentum and the direction opposite to S s i-momentum in the D*>** 
rest frame. To extract the D*, Dq and the D\ helicity distributions we perform a combined 
fit of the M(Dtt) spectra for combinations from both B + and B° in bins of helicity 
angle. The fit procedure is identical to that used for the B(B — > D*'**£u) calculation. 
The results corrected for the efficiency are plotted in Fig. [31 D\ distributions for D* v and Dq 
hypothesis coincide within errors, so that only that for the Dq +D% case is shown in Fig. [3]c. 
The .Dq helicity distribution is consistent with the J = hypothesis (x 2 /ndf = 6.0/4, 
where ndf is the number of degrees of freedom). The D% helicity distribution is fitted 
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TABLE II: Results of the D^tt+ pair invariant mass study. i3(mode) = B(B -» D**£u)xB{D** -> 
L»(*)vr + ). The first error is statistical and the second is systematic. 



Mode 


Yield 


£(mode),% 


Signif. 


B+ 


u 


102 ± 19 


0.24 ±0.04 ±0.06 


5.4 


B+ 


-> Dft+v 


94 ± 13 


0.22 ±0.03 ±0.04 


8.0 


B°- 


-4 D*-£+u 
u 


61 ±22 


0.20 ±0.07 ±0.05 
< 0.4 @ 90% C.L. 


2.6 


B°- 




68 ± 13 


0.22 ±0.04 ±0.04 


5.5 


B+ 


-> s^+i/ 


— 5± 11 


< 0.07 @ 90% C.L. 




B+ 




81 ± 13 


0.42 ±0.07 ±0.07 


6.7 


B+ 




35 ±11 


0.18 ±0.06 ±0.03 


3.2 


B° 




4±8 


< 0.5 @ 90% C.L. 




B° 


-» L>r 


20 ± 7 


0.54 ±0.19 ±0.09 
< 0.9 @ 90% C.L. 


2.9 


B°- 


-4 D^+z/ 


1 ±6 


< 0.3 @ 90% C.L. 





with the function c^o I ^2° 1 2 + 4af | V^ 1 1 2 + 4a| 1 1 2 , where the Y- are spherical harmonics and 
al + 4a 2 + 4al = 1. The fit yields a 2 = 0.74 ± 0.10, a 2 = 0.04 ± 0.02 and a\ = 0.02 ± 0.02; 
the fit quality is x 2 /ndf = 2.0/3. The fit is consistent with the assumed quantum numbers 
and demonstrates that the D 2 from semileptonic decay is dominantly in the s z = spin 
projection. Helicity distributions, predicted by theory, are shown as dashed lines. For 
evaluating the D* + D 2 hypothesis, the obtained D* helicity distribution (Fig. [3] b) is fitted 
with the function b 2 \Y®\ 2 + b\\Y^\ 2 . This fit yields b 2 = 0.15 ± 0.09, b 2 = 0.85 ± 0.09 
(x 2 /ndf = 18.8/4) in poor agreement with expectations from theory, shown as a dashed 
line. 

We also study the dependence of the B — > D** transition on q 2 or, equivalently, on 
the conventional HQET variable w, which is the dot-product of B and D** four-velocities: 
w = vb ■ Vd**- The ^-dependence is obtained from fits of Dtt invariant mass in bins of w. 
The results are presented in Fig. |U As with the helicity study the D 2 distribution is shown 
only for the Dq + D * hy pothesis in Fig. H] c. The w distribution is fitted according to the 
model given in Ref. [17|. In HQET, the matrix elements between the B and D states to 
leading order in Aqcd/^q are expressed in terms of three universal Isgur-Wise functions 
f(tu), t 1/2 (w) and t 3/2 (w) for (D, D*), (D%,D[) and (D U D%) doublets, respectively [[17]. We 
assume a "pole" form for £(w): £ = (2/(1 + w)) 2p2 and a linear form for Tj(to) functions: 



Ti(w) = Tj(l)[l +fl(w — 1)], and the following relation: f[, 2 = t' z , 2 + 0.5 [18j. A simultaneous 
fit to the ^-distributions for and D 2 gives fy 2 = —1.8 ±0.3. Using the measured branch- 
ing ratios of B — > Dq 2 £u, we also calculate t 3/ / 2 (1) = 0.75 and Ti/ 2 {1) = 1.28. All parameters 
are in agreement with expectations except for 71/2(1), which is larger than predicted due to 
the large value of B(B — > D^lv). 

The systematic error in the calculation of branching fractions due to B t3ug efficiency un- 
certainty cancelled out since normalization modes were used. The largest contribution to 
the systematic error is from uncertainty in D** parameters. We perform a D^*'n mass study 
with these parameters allowed to float inside their errors to get it. The B s \ reconstruction 
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efficiency dependence due to the decay model was studied using two different signal MC 
samples generated with the ISGW2 fj| and Goity-Roberts [l6[ models. To estimate the 
systematic uncertainty in background subtraction we used two different sets of sidebands 
with appropriate normalizations. To estimate interference effects we perform MC study with 
different angle efficiency dependencies. A summary of the systematic error contributions is 
presented in Table Hm In total we obtain a 14% error for the B — > D^niu measurement, 
a 16% error for the narrow D** contribution and 25% for the measurement of the wide D** 
contribution. 



TABLE III: Sources of systematic error. 



Source of error 


Systematic error, % 


Reconstruction 


2 


Model efficiency 


7 


Background subtraction 


6 


Misidentified hadrons 


1 


Interference 


7 


Normalization branching ratios 


10 


D** parameters 


6-20 


Total 


14-25 



In conclusion, we report measurements of the branching fractions for B — > D^niu decays. 
These measurements supersede our previous results (8|. We also performed an analysis of 
the final state D^'ir hadronic system and obtained branching ratios for the B — > D**£u 
components. Semileptonic decay to D\ meson is observed and measured for the first time. 
Helicity and w distributions are studied for this decay. We observe a broad enhancement 
in the Dir mass distribution consistent with wide Dq production. The branching ratio of 
the decay to B — > D^lv is found to be large, in contrast with theoretical predictions |4j. 
However there is no indication of a broad D' x in the B — > D*ir£i> channel, which should be of 
the same order. The combined likelihood of fits to the D-ir mass, helicity and w distributions 
for Dq + D\ hypothesis is higher than that for the D* v + D\ combination by 2.8a. 

However, the present data sample cannot exclude the interpretation of this enhancement 
as a D* tail. 
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FIG. 1: Ml spectra before (1) and after (2) background subtraction for: a) B + — > D ir + l + v, b) 
B+ D*-n + £ + v, c) B° -► D°7t~£ + v, d) B° -► D* ^-l + v. The curves are the fits, which are 
described in the text. 
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FIG. 2: Hadronic invariant mass distributions for: a) B + — > Z) 7r + ^ + z/, 6j i? + — > D* 7r + £ + ^, 

7T ^ + z/, (fj £?° — > Z?*°7r .£ + za Insets show the distributions before background subtraction 
in the region around the narrow D**'s. The background is shown as the hatched histogram. The 
curves are the fits, which are described in the text. 
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FIG. 3: Helicity distributions for a) Dq, b) D*, c) D\. The curves represent the fits, described in 
the text. 




FIG. 4: w distributions for a) Dq, b) D*, c) D\. The curves are the fits, which are described in 
the text. 
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